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ABSTRACT 



AGN feedback now appears as an attractive mechanism to resolve some of the outstanding problems with the "standard" 
cosmological models, in particular those related to massive galaxies. At low redshift, evidence is growing that gas cooling 
and star formation may be efficiently suppressed by mechanical energy input from radio sources. To directly constrain 
how this may influence the formation of massive galaxies near the peak in the redshift distribution of powerful quasars, 
z~ 2, we present an analysis of the emission-line kinematics of 3 powerful radio galaxies at z~ 2 — 3 (HzRGs) based 
on rest-frame optical integral-field spectroscopy obtained with SINFONI on the VLT. The host galaxies of powerful 
radio-loud AGN are among the most massive galaxies, and thus AGN feedback may have a particularly clear signature 
in these galaxies. 

We find evidence for bipolar outflows in all HzRGs, with kinetic energies that are equivalent to 0.2% of the rest-mass 
of the supermassive black hole. Observed total velocity offsets in the outflows are ~ 800 - 1000 km s"^ between the 
blueshifted and redshifted line emission, and FWHMs ~ 1000 km s~^ suggest strong turbulence. Line ratios allow to 
measure electron temperatures, ~ 10* K from [OIII]AAA4363, 4959, 5007 at z~ 2, electron densities (~ 500 cm~^) and 
extinction (Ay ~ 1 — 4 mag). Ionized gas masses estimated from the Ha luminosity are of order 10^" M©, similar to 
the molecular gas content of HzRGs, underlining that these outflows may indicate a significant phase in the evolution 
of the host galaxy. The total energy release of ~ lO'''^ erg during a dynamical time of ~ 10^ yrs corresponds to about 
the binding energy of a massive galaxy, similar to the prescriptions adopted in galaxy evolution models. Geometry, 
timescales and energy injection rates of order 10% of the kinetic energy fiux of the jet suggest that the outfiows are 
most likely driven by the radio source. The global energy density release of ~ lO^'^ erg s~^ Mpc~^ may also influence 
the subsequent evolution of the HzRG by enhancing the entropy and pressure in the surrounding halo and facilitating 
ram-pressure stripping of gas in satellite galaxies that may contribute to the subsequent mass assembly of the HzRG 
through low-dissipation "dry" mergers. 

Key words. Keywords should be given 
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^ 1. Introduction 

AGN feedback has now become a crucial process in most 
models of galaxy formation and evolution. Current mod- 
els postulate that massive galaxies experienced an early 
epoch of vigorous star formation that was terminated by a 
nearly instanta neous, power ful bl ow-out phase triggered by 
Siik fc Reeslll998l: IScannapieco &: Ohll2004t 



the AGN (e.g. 
Springel et al.1 



20051: iDi Matteo et al.ri2005l: ICroton et all 



2006tlBower et al.ll2006l:lHopkins et al.H2006[ ). Such a meek 



anism would heat and remove most of the ambient gas of 
the galaxy, and thus reconcile the discrepancy between the 
kierarchical standard model of galaxy formation and obser- 
vational constraints 

Bolometric luminosities of powerful AGN at kigk red- 
skift correspond to an energy output of 10^°"^^ erg during 
an AGN lifetime of 10^~^ yrs. Altkougk tkis is equal, or 



* Based on observations collected at the European Southern 
Observatory, Very Large Telescope Array, Cerro Paranal, Chile 
( 076.A-0684(A)) 



even exceeds the binding energy of a massive galaxy, tkis 
energy will only kave an impact on tke evolution of tke kost 
galaxy, if it is efficiently transferred to tke ambient gas witk 
a coupling efficiency of at least a few percent. 

Observationally, evidence is growing tkat AGN feed- 
back may be related mostly to radio-loud AGN. Perkaps 
tke most spectacular examples at low redskift are tke gi- 
ant cavities in tke X-ray kalos of massive galaxy clusters, 
wkick appear to inject few x 10^" erg into tke surrounding 
gas witk in a few xlO^ y r s, sufficient to suppress gas cooling 
(Bockri nger et all 119931: iBirzan et all 120 04; Raffert v^eTall 
i2006c ,McNamara fc Nulsenll2007[ ). iBest et al.l (|2005t l2006f) 
find more subtle, yet general evidence tkat radio sources 
may balance gas cooling in ~ 2000 radio-loud early-type 
galaxies taken from tke Sloan Digital Sky Survey and tke 
NVSS and FIRST radio surveys. Interestingly, tke frac- 
tion of radio loud galaxies seems to be a strong function 
of mass, suggesting tkat suck radio-driven feedback will 
act predominantly in tke most massive galaxies, in broad 
agreement witk tke predictions of galaxy evolution mod- 
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els. This is also supported by the observation that pow- 
erful radio galaxies are among the most massive galaxies 
at all red shifts (De Breuck et al. 2001; Willott et al. 2003; 
lRocca-\^lmerange et al.ll2004f ). Since AGN feedback is ex- 
pected to act predominantly in massive galaxies, we may 
expect that its influence will be most clearly expressed in 
particularly massive galaxies. 

However, since massive early- type galaxies appear to 
have completed mos t of their growth at high redshift (e.g., 
iRudnick et al.l I2003L and references therein), directly in- 
vestigating the impact of jet-driven AGN feedback dur- 
ing the formation of massive galaxies requires direct ob- 
servations of radio-loud galaxies at z> 1. Interestingly, 
longslit spectroscopy of powerful radio galaxies at z> 1 
(and in a few cases at lower redshifts) revealed strongly 
distorted emission line kinematics, with velocities often 
exceeding ^ 1000 km s~ ^ along the axis of th e radio 



jet ( e g., [Tad huntei J [199]]: iMcCart hv et al.l [Till 



19981 iBaum fc McC arthvl 120001: [yjUar-Ma rtm et 8,1 
Inskip et alf l2002i : iClark et all \l998i iBest et al 



■ 19991 : 
19971) . 



Similarly, [OIIIJA5007 narrow-band imaging of 4C19.71 at 
z~ 3.6 reveale d a giant emission line region extending 
over ~ 74 kpc (jArmus et al.l [19981 ). This emission line re- 
gion, corresponding to an ionized gas mass of 10^~^ Mq 
aligns rough ly with the axis of the radio jet, which led 
lArmus et al.l to suspec t that the jet and the neb ulosity may 
be physically related. iHeckman et al.l (|1991bl lal) also find 
this for radio-loud quasars. Moreover, Lya spectroscopy 
and line imaging revealed luminous emission line halos with 
sizes that appear related to the size of the radio jet, some- 
times with kinematically more quies cent gas beyond (e.g. 
IVillar-Martm et al.ll2003L I2006L I2007D . 

However, the complexity of the velocity fields and rel- 
atively large spatial extent of the emission line regions 
made it difficult to robustly infer the global characteris- 
tics of these nebulae. Integral-field spectrographs now open 
a new avenue to study the underlying dynamical mecha- 
nisms with unprecedented robustness and across the full 
two-dimensional surface of the emission line nebulae. Over 
recent years, our understanding of galaxy evolution in the 
early universe has significantly advanced due to the dis- 
cover y of large numbers of galaxies at similar redshifts 
(;e.g.. lSteidel et al.lll996tlSmail et al.ll2002D . and during this 
time, modeling of galaxy and structure formation has im- 
proved. This allows us to reassess the emission line kine- 
matics of powerful high-redshift radio galaxies (HzRGs). 
Near-infrared studies are particularly well suited to study 
the inner regions of HzRG extended emission line regions, 
allowing for relatively high spatial resolution of ~ 0.5" even 
with seeing-limited observations, and opening a window to 
the rest-frame optical line emission where extinction is less 
se vere than in the r e st-fram e UV. 

iNesvadba et al.l (|2006bl ) find strong evidence for en- 
ergetic outflows of ionized gas from rest-frame optical 
integral-field spectroscopy of the powerful radio galaxy 
MRC1138-262 at z~ 2.2. With a spatial extent of 30 kpc, 
relative velocities and line widths of FWHM~ 1000 km 
across the source, corresponding to kinetic energies 
of few X 10™ erg, this gas may be experiencing a feedback 
episode as dramatic as necessary to unbind a significant 
fraction of the ambient gas from the halo of a massive 
galaxy. The kinetic energy of the outflow corresponds to 
a few percent of the jet kinetic energy, and the dynam- 
ical timescale of a few xlO^ yrs appears roughly similar 



20031) : iDe Breuck et all fe.g.. l2000D:lHumDhrev et al. 
2008f) ; IBaum fc McCarthvl fe.g.. l2000f) . we find that 



to the age of the radio source. Moreover, the turbulent, 
luminous gas extends to the size of the radio jets. In con- 
trast, very similar observations do not show evidence for 
such outflows in compact r adio galaxies at similar red- 
shifts (jNesvadba et al.ll2007aD . which are likely the younger 
analogs to extended HzRGs. This agrees with expectations. 
If the entrainment rates are similar, then only a few percent 
of the ISM will already be affected by the radio source for 
ages less than '^10^ yrs. 

Currently, the host galaxies of powerful, radio-loud 
AGN (HzRGs and quasars) represent the only galaxy popu- 
lation at z> 2 with such extreme kinematics. To address the 
question whether MRC 11 38-262 is a "one-off", or whether 
AGN feedback is a c ommon phenomenon in HzRGs , we con- 
tinue the analysis of' Nesvadba et al.l (|2006bl . [2007a[ ). adding 
another three HzRGs with near-infrared int egral-field data 
at 2 — 3 . Comparing with the s amples of Iwamuro et al 
(e-g-. 
(e-g-. 

their rest-frame UV-optical properties are within the typi- 
cal range of powerful HzRGs. These galaxies have extended 
radio morphologies, and relatively similar radio power to 
MRC1138-262, Viaghz ~ 10^^ W Rz'K Thus, they are 
among the most powerful radio galaxies observed at all 
redshifts. Given the observational difficulties in studying 
high-redshift galaxies, this seems the most promising way 
of identifying the fingerprints of powerful AGN-driven feed- 
back and underlying physical processes directly from the 
gas kinematics of powerful AGN host galaxies at z^ 2. 

However, th ese radio powe r s are not exceptional in the 
early universe. IWillott et el\ ()200lD studied the redshift 
evolution of the radio luminosity function and find a popu- 
lation of galaxies with radio powers Vi^imhz ^ io^^'^~^^'^ 
W Hz^^ at 151 MHz, which shows a strong peak near z~ 2 
and rapidly declines towards lower redshifts. Finally, the 
radio emissivity is likely a strong function of environment 
and evolutionary s tage of the synchrotron plasma (e.g., 
iKaiser et al1ll997al [bl). causing changes in the radio lumi- 
nosity of more than an order of magnitude. Hence, by se- 
lecting particularly powerful sources, we may bias our sam- 
ple more in terms of the evolutionary stage or environment 
than the kinetic luminosity of the radio source. 

Specifically, we observed MRC03 16-257 at z= 3.13, 
MRC0406-244 at z= 2.42, and TXS0828-I-192 at z= 2.57. 
Entrained ionized gas masses and kinetic energies seem very 
similar to MRC1138-262 in all three cases. This leads us to 
propose that AGN feedback may be a common phenomenon 
in powerful radio galaxies during the "Quasar Era" . 



2. Observations and Data Reduction 

Observations in the H and K bands were ca rried out with 
the integral-field spectrograph SINFONI ([Bonnet et al.l 
1200 4) on the VLT in December 2005 under good conditions. 

SINFONI is a medium-resolution, image-slicing 
integral- field spectrograph, with 8"x 8" field of view 
at a 0.25" X 0.25" pixel scale, and spectral resolution 
of i? ~ 3000, and R - 4000 in the H, and K bands, 
respectively. Individual exposure times were 450 s in H and 
600 s in K, respectively. Total exposure times are given for 
each band and each so urce in Table [TJ 

We use the IRAF (|Ibdvl [l99^ standard tools for the 
reduction of longslit-spectra, modified to meet the spe- 
cial requirements of integral-field spectroscopy, and comple- 
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mented by a dedicated set of IDL routines. Data are dark- 
frame subtracted and flat-fielded. The position of each slit- 
let is measured from a set of standard SINFONI cahbration 
data, measuring the position of an artificial point source. 
Rectification along the spectral axis and wavelength cali- 
bration are done before night sky subtraction to account for 
some spectral flexure between the frames. Curvature is mea- 
sured and removed using an arc lamp, before shifting the 
spectra to an absolute (vacuum) wavelength scale with ref- 
erence to the OH lines in the data. To account for variations 
in the night sky emission, we normalize the sky frame to the 
average of the object frame separately for each wavelength 
before sky subtraction, masking bright foreground objects, 
and correcting for residuals of the background subtraction 
and uncertainties in the flux calibration by subsequently 
subtracting the (empty sky) background separately from 
each wavelength plane. In the H band, where the suppres- 
sion of night sky lines is particularly difficult, we addition- 
ally clipped all strongly deviant pixels at wavelengths near 
prominent night sky lines. 

The three-dimensional data are then reconstructed and 
spatially aligned using the telescope offsets as recorded in 
the header within the same observing block (a sequence of 
typically 5 to 10 individual exposures with a total length 
of about an hour), and by cross-correlating the line im- 
ages from the combined data in each observing block, to 
eliminate relative offsets between different observing block. 
Telluric correction is applied to each individual cube be- 
fore the cube combination. Flux scales are obtained from 
standard star observations. From the light profile of the 
standard star, we measure the FWHM spatial resolution in 
the combined cube, which is typically ^ 0.6". 

2.1. Complementary Data Sets and Alignment 

We complement our SI NFONI data throug h deep, A-array 
1.4 GHz radio maps of ICarilli et al.l (|1997 ) with a spatial 
resolution of ^ 0.4". These data were kindly provided to 
us by C. Carilli. Studying the interplay between radio jets 
and the ISM of high-redshift galaxies requires a relative 
alignment that is accurate to < 0.5", which is better than 
the accuracy of the absolute astrometry reached with the 
VLT. We therefore assume that the radio core falls onto the 
center of the continuum emission for all sources. Note that 
this coincides with the position of the largest velocity gradi- 
ent in the velocity maps of TXS0828-f 193 and MRC0406- 
244 to within less than a few tenths of a second of arc. 
Remaining small offsets are likely due to projection effects. 
Since we did not detect MRC0316-257 in the continuum, 
we align the radio core with the steepest gradient in the 
velocity map, implying that the basic physical mechanism 
of the outflow is similar in the 3 sources - which is justi- 
fied given their ov erall similarity. We also note that, unlike 
in MRC1138-262 (|Nesvadba et al.ll2006bD . we do not find a 
nuclear broad-line region related to the AGN in any of the 
three sources. 



3. Continuum and Line Emission Morphologies and 
Integrated Spectra 

We constructed continuum-free line images by summing 
over the spectral regions of the emission lines in the in- 
tegrated spectra, and removing the underlying continuum 




2,00 2.02 2.04 2.06 2.08 
obs. wovelength [f^m] 



Fig. 1. Integrated K band spectrum of MRC0316-257 
at z=3.13, extracted from all pixels in which the 
[OIIIJA5007 line emission exceeded 3 a, and covering the 
[OIII]AA4959,5007 doublet and H/3. 

emission. Rather than adopting a fixed aperture, we se- 
lected all pixels where the [OIIIJA5007 emission line was de- 
tected at a signal-to- noise ratio greater than 3. This helps to 
maximize the signal-to-noise ratio, which may otherwise be 
affected by night-sky residuals, whose strength often varies 
across the detector. To extract line- free continuum images, 
we fitted the spectrum of each spatial pixel where it is un- 
affected by bright emission lines or strong night sky line 
residuals. We interpolated over the contaminated parts be- 
fore collapsing the spectra. 

3.1. MRC0316-257 

The integrated spectrum of MRC0316-257 is shown in 
Fig.[U The [OIII]AA4959^5007 doublet and H/3 are detected 
at redshifts of z=3. 127^1 and z=3.1292 for [OIII]A5007, re- 
spectively, and have line widths of FWHM=1465 km s~^ 
and FWHM=900 km s^^ for [OIII]A5007 and H/3, respec- 
tivel}0. See Table [2] for a summary of the emission line 
properties and related uncertainties. 

We show the [OIIIJA5007 emission line morphology in 
Fig. [51 Line emission is detected over an area of w 2 . 8" x 1 . 2" 
(21 kpc x9 kpc), and is elongated along the axis of the radio 
jet (the 1.4 GHz VLA map is shown as contours in Fig. [5]). 
H/3 and continuum emission are very faint, so that their 
morphology could not be measured. 

3.2. MRC0406-244 

MRC0406-244 at z=2.4^has a splendid hour-glass shaped, 
rest-frame UV and optical morphology which has previ- 

^ Adopting the flat f^A ~ 0.7 cosmology with Hq = 70 km 
Mpc~^ we use a luminosity distance Dl = 26.7 Gpc and 
angular size distance Da = 1.5 Gpc. The size scale is 7.6 kpc/". 
The age of the universe at z~ 3.13 for this cosmological model 
is 2.0 Gyr. 

^ Unless stated explicitly, all kinematic parameters are in the 
rest-frame. 

^ The luminosity distance is Dl = 19.8 Gpc and angular size 
distance Da = 1-7 Gpc. The size scale is 8.1 kpc/". The age of 
the universe for this redshift and cosmological model is 2.6 Gyr. 



4 



Nesvadba et al.: Jet-driven AGN winds in Radio Galaxies in the "Quasar Era" 



8000 



6000 



i 4000 



2000 



[QUI] 




"'J'J'KAyU'W' 



1.55 1,60 1.65 1.70 1-75 

obs. WQveJength [^m] 



4000 



3000 r 



1000 r 




lOOO 



Ha+[NII] 

II 



2.J-> 2.26 2.26 2.30 2.J2 



[SII]67IS,S73t 
II 




2.00 3.10 2.20 2.30 

obs. wavelength [p-m] 

Fig. 2. Integrated H (top) and K band spectrum (bottom) 
of MRC0406-244 at z=2.42, extracted from all pixels in 
which the [OIIIJA5007 line emission exceeds 3 a. The in- 
set in the upper panel shows the detected, relatively faint 
H/3, and [HeII]A4685 lines, and the expected wavelength of 
[HeI]A4471, which is not detected. The inset in the right 
panel shows (from top to bottom): (1) the blended Ha and 
[NII]AA6548,6583 hues, and the [SII]AA6716,6731 doublet 
(2) the best Gaussian fit to these lines (see text) (3) fit 
residuals. The fitted lines and residuals are shifted along 
the ordinate by arbitrary amouts. 



ously been studied with NICMOS a nd WFPC imaging on - 
board the Hubble Space Telescope (iPentericci et al. 200lh. 
as we l l as w it h ground-based ob servations (|McCarthv et al.l 
119911 119921 : iRush et alj I1997D in the optical and near- 
infrared. As pointed out bv iRush et all (|l997t ). the mor- 
phology is a strong function of waveband. In particular, 
the J-band morphology is much m ore compact than, e.g., 
the morphology in the Ks or g-band. TRush et al] (jl997[ ) sus- 
pected that this is most likely due to contamination with 
strong emission lines (Ha in the Ks band, Lya in the g 
band). 

Our SINFONI data illustrate that this is indeed the 
case. In Fig. IHlwe show the [OIIIJA5007 emission line mor- 
phology. Contours indicate the line-free K-band continuum. 
Line emission is clearly extended over 3.7" x 1.1" (30.1 
kpcx 9.3 kpc) and is reminiscent of two edge-brightened 
bubbles with severa l brighter knots super imposed. HST 
NICMOS imaging of IPentericci et all ()200lD in the F160W 



filter (which covers [OIIIJA5007 and is dominated by line 
emission) reveals a more complex morphology, with several 
bright knots and fainter filaments inbetween. 

The line-free continuum morphology extracted from the 
SINFONI data cube shows that continuum emission is faint 
and relatively compact, but spatially resolved (contours in 
Fig. [5]). We detect the continuum over an area of approx- 
imately 1.3"x0.6" (10.6 kpcx5 kpc), corresponding to an 
integrated magnitude of K= 17.7 mag in a 3" aperture cen- 
tered on the continuum peak (compared to K= 17.5 mag if 
not correcting for the line contamination). In the H band, 
we measure H= 18.9 mag and H= 19.3 mag for the uncor- 
rected and corrected magnitude, respectively. We also esti- 
mate an azimuthally-averaged half-light radius, which cor- 
responds to the area of brightest continuum emission, that 
contains half the total continuum fiux above a 3a threshold. 
We find a half-light radius of ri/2 ~ 0.6" 5 kpc). 

The integrated spectra of MRC0406-244 in the H and 
K bands are shown in Fig. [51 Rest-frame optical emission 
lines are very bright, in particular the [OIII]AA4959,5007 
doublet and H/3 in the H band, and the blended Ha and 
[NIIJA6583 emission lines in the K band. In the K band, 
we also detect the [SII]A6716,6731 doublet and [OI]A6300. 
Moreover, we marginally detect the faint HeIIA4686 in the 
H band spectrum. The inset in the left panel of Fig.[2]shows 
the fainter lines in the H-band. The profiles of some lines are 
aflfected by residuals of night-sky lines (e.g., on the red wing 
of H/3). See TableOfor the observed wavelengths, redshifts, 
line widths and fluxes of the detected emission lines in the 
integrated spectrum. 

3.3. TXS0828+193 

TXS0828+193 at z=2.5'fl is the galaxy with the most ex- 
tended radio lobes in our sample (103 kpc for our cosmol- 
ogy). F675W WFPC2 imaging indicates a clumpy, irregular 
and asymmetric morph ology at rest-frame UV wavelengths 
(Pentericci e t al.ll200lD . 

TXS0828-I-193 is a particularly luminous line emitter. In 
Fig. [H we show the integrated spectrum of TXS0828-hl93 in 
the H and K bands. [OIII]AA4959,5007, H/3, H7 are detected 
in the H band. [OIIIJA4363 is faint, but detected with a sig- 
nificance of 3.6 a. In the K band we observe [OIJA6300, Ha 
and [NII]A6583. The [SII]AA6716,6731 doublet falls at ob- 
served wavelengths A>2.4, which are dominated by telluric 
absorption. The emission line properties of the integrated 
spectrum of TXS0828+193 arc summarized in Table |4l 

We show the [OIIIJA5007 emission line morphology of 
TXS0828-f 193 in Fig.[TOl The bright extended emission line 
regions are significantly more extended than the continuum 
(shown as contours in the left panel of Fig. [TU)) . covering 



an area of 3.5"x 1.2" ( 28 kpcx9.6 kpc). The line image is 
dominated by two bright regions, an unresolved peak to the 
south, and an extended, more diffuse emission line region to 
th e north-east. Compariso n with the WFPC2 F675W image 
of Pcntericci et aTl (|2001f ) (their Fig. 7) indicates that the 
[OIIIJA5007 emission line morphology and the rest-frame 
UV morphology are very similar. 

The continuum is faint in the H and K band (approx- 
imately corresponding to the rest-frame V and R bands) 



The luminosity distance is Dl = 21.1 Gpc and angular size 
distance Da = 1-7 Gpc. The size scale is 8.0 kpc/". The age of 
the universe for this redshift and cosmological model is 2.5 Gyr. 
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Fig. 3. Integrated H (top) and K band spectrum (bottom) 
of TXS0828+193 at z=2.57, extracted from all pixels in 
which the [OIIIJA5007 line emission exceeds 3 a. The inset 
in the upper panel shows the detected, relatively faint H/3, 
H7, and [HeII]A4685 lines, and the expected wavelength of 
[HeI]A4471, which is not detected. [OIII]AA4959,5007 have 
complex profiles, where the different components arise from 
the blue and redshifted bubble, respectively. [OIIIJA6300 
is faint, but detected with a significance of 5cr. In the K 
band, we detect Ha and [NII]AA6548,6583, and the fainter 
[OI]A6300 line. The [SII]AA6716,6731 doublet is at observed 
wavelengths > 2.4/zm, which are dominated by telluric ab- 
sorption. 



and significantly more compact than the line emission. We 
show the line-free, combined H and K band continuum 
morphology as contours in Fig. 1101 Continuum emission 
is detected over an area of ^ 1.4"x0.8", corresponding 
to about 11.3 kpcx6.4 kpc. Estimating an azimuthally- 
averaged half-hght radius similar to ^3.21 we find 0.47" (3.8 
kpc). region. We estimate integrated magnitudes in the H 
and K band using 3" apertures, K=19.4 mag, H=19.7 mag 
(compared to K=18.5 mag, and 11=19.1 mag, without cor- 
recting for line contamination). 



4. Characteristics of the extended emission line 
regions 

4.1. Kinematics 

We extracted spectra from 3x3 pixel square apertures 
across the sources (0.375" x 0.375") and fitted Gaussian line 
profiles to each spectrum. Relative velocities were estimated 
from the measured wavelengths of the line cores at each po- 
sition (core of the strongest component in case of complex 
line profiles). 

The large line widths observed in the extended emission 
line regions of HzRGs complicate the analysis of the spectra 
due to the blending of emission lines, in particular of Ha 
and [NII]A6548,6583. In order to obtain a consistent fit of 
all lines within each spectrum, we therefore estimated the 
redshift and line width from the [OIII]AA4959,5007 doublet 
(whose components are luminous and can be isolated rela- 
tively easily) and imposed the same redshift and line width 
on all other emission lines. As a result, the line fluxes are 
the only free parameter in the fits to all other lines except 

[oiii]. 

However, such an approach implies that all emission 
lines trace the same gas kinematics. We tested this explic- 
itly by comparing velocities and line widths obtained from 
[OIII]AA4959,5007 in TXS0828-)-193 with those measured 
from H/3, in areas where H/3 is detected at a significant 
level (> 3cr), and find good agreement between the mea- 
sured velocities and line widths. 

For much of the extended emission, the residuals sug- 
gest that fitting a single Gaussian component is ade- 
quate. However, in the central regions in particular, we 
do observe more complicated profiles and line splitting 
(see, e.g.. Fig. [7]). In these cases, we fitted 2—3 compo- 
nents with redshifts and line widths determined from the 
[OIII]AA4959,5007 doublet, and leaving the flux of each line 
and each component as a free parameter. 

We present the maps of rest-frame relative velocities and 
line widths (full width at half maximum) of the 3 galaxies 
in Fig. [5] to [TO] . Overall, the velocity fields are remarkably 
similar in all targets. In each galaxy we identify two regions 
with relatively homogeneous internal kinematics, which ex- 
tend from near the radio core, and projected velocities rela- 
tive to each other of ~ 700—1000 km s~^. Velocities change 
abruptly in the central region of the velocity field, within 
about a spatial resolution element of our data sets, corre- 
sponding to ^ 0.6" or < 5 kpc. The velocity gradients for 
the individual galaxies are given in Table [HI 

The velocity maps show substructure within each re- 
gion, although the velocity gradients are significantly 
smaller than the large velocity shift near the center of tlie 
galaxies. Velocity gradients amount to ^ 100—200 km s~^. 
We emphasize that the three targets have similar velocity 
fields, in spite of different morphologies. 

Figures [5] to [TO] also show maps of the FWHM line 
widths. Overall, line widths are FWHM-^ 500 - 1200 km 
s^^ in the three galaxies. Note that much of this range 
of FWHMs is due to variations in the line widths within 
individual areas, and not so much due to large variation 
between galaxies in our sample. 

Apparent peaks in the central regions with FWHM> 
1400 km s~^ are likely artefacts due to the large velocity 
gradients and partial overlap of line emission from the two 
sides within one seeing disk. We also observe line splitting 
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(see Fig. [7]) in this area, and the observed line widths are 
consistent with what would be expected from blended lines 
with the observed line widths and velocity gradient in and 
between the two regions. In other regions of the galaxies, 
the emission line widths likely reflect a large intrinsic turbu- 
lent motion, since the widths exceed the velocity gradients 
within each bubble (at least on kpc scales). 

4.2. Extinction 

MRC0406-244 and TXS0828-hl93 are bright Ha and H/3 
emitters and allow to construct maps of the H/J/Ha line 
ratios and thus extinction. For MRC0316-257 this is not 
possible, due to the larger redshift, z=:3.13, at which Ha 
cannot be observed from the ground, and because H/3 is 
relatively faint. 

We measured the H/3/Ha line ratios from 0.375" x 
0.375" box apertures (3x3 pixels). To derive extinctions, 
we use a galactic extinction law and an intrinsic (unred- 
dened) Balmer decrement Ha/H/3 = 2.9. (See i )4.1l for a 
description of how we correct for blending between Ha and 
the [NH]AA6548,6583 doublet.) 

We show the extinction maps of MRC0406-244 and 
TXSG828-f 193 in the top and bottom panel of Fig. H re- 
spectively. Contours indicate the [OHIJA5007 emission line 
morphology to ease orientation. Extinctions at the wave- 
length of H/?, A{Hf3) vary between A{H(3) ^ 1 — 4.5 mag 
in both galaxies. (We give A(H/3)=A(V)+0.14 mag, be- 
cause we measured the extinction from the emission lines.) 
The most heavily extincted regions are near the center of 
each source, where A(Hf3)^ 4 — 4.5 ; the spatial extent 
of these areas corresponds approximately to the size of 
the continuum emission. Extinction in the extended emis- 
sion line regions is less severe, with A{H(3)^'^^^ 1 — 2 
mag in MRC0406-244, and A{HP)°lf ~ 2 - 3 mag in 
TXSG828+193. 

Actively star-forming galaxies may have strong stellar 
Balmer absorption lines, in particular for the higher-order 
Balmer lines. This may result in underestimating the H/3 
emission line flux relative to Ha, and thus in overestimat- 
ing the extinction. Since the continuum in our galaxies is 
faint we cannot measure the effect of underlying stellar ab- 
sorption directly. However, the large Ha and H/3 equivalent 
widths, the faintness of the continuum, and the significantly 
larger spatial extent of the line emission indicate that stel- 
lar absorption will have a negligible effect. 



4.3. Electron Densities 

The [SH]AA6716,6731 doublet in TXS0828+193 is relatively 
faint and falls into a wavelength range that is dominated 
by telluric absorption bands; however, in MRC0406-244, 
the doublet is clearly detected in the K band at signal-to- 
noise ratios « 10. (see Fig. [5] for the integrated K band 
spectrum of MRC0406-244). Due to the large fine widths, 
the doublet is partially blended. To minimize the blend- 
ing, we extracted the integrated spectrum of both bubbles 
individually, avoiding the central region where lines are par- 
ticularly broad due to the overlap of the two bubbles. This 
is justified by the overall very uniform line ratios of the 
[SII]AA6716,6731 doublet extracted from smaller apertures. 
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Fig. 4. A(H/3) extinction maps for MRC0406-244 (top) 
and TXSG828-hl93 (bottom). Color bars show the extinc- 
tions (in magnitudes) estimated from the observed H/3/Ha 
emission line ratios, assuming a galactic extinction law. 
Contours show the emission line morphology to ease ori- 
entation. 



We show the K band spectra of the two bubbles in 
Fig. [S] The upper panels show the measured spectra, fits 
and fit residuals are shown below. The profiles of all lines 
in the K-band are relatively narrow and regular. We used 
again [OHI]A4959,5007 in the H band to constrain the red- 
shift and line width, so that the line fluxes are the only free 
parameters. 

[SII] line ratios, R[sh] = F(6716)/F(6731), in the two 
bubbles are similar within the uncertainties; we measure 



R|'406|. ^ 1.0 ±0.15 in the blueshifted bubble. 



and R«|??f 
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Fig. 5. Integrated spectra of the two bubbles in MRC0406- 
244, the redshifted southern bubble is shown on top, the 
blueshifted northern bubble is shown below. Upper pan- 
els show the spectra, Gaussian fits to the emission lines 
and fit residuals are shown in the mid and lower panels, 
respectively. Spectra are smoothed by 5 pixels spectrally, 
corresponding to ^ 1/5 of the total line width. The mea- 
sured [SII]AA6716,6731 line ratios are consistent with elec- 
tron densities of a few 100 cm~^, with a best-fit value of 
500 cm~^ (see text). Redshifts and line widths were mea- 
sured from the [OIII]AA4959,5007 hne doublet, so that the 
only free parameters to all other lines, including [SII], are 
the line fluxes. 
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1.0 ± 0.10 in the redshifted bubble. For a temperature of 
T= 10*^ K (Sm, this corresponds to a range in electron 
density ofn°*°^'' = 300-1045 cm'^ and nO^osr = 370-830 
cm~^ in the blueshifted and redshifted bubble, respectively, 
with a best value of about 500 cm~'^ in either case, which 
is the value we will adopt in the following. 

Only one z~ 2 HzRG has an electron density es- 
timate given in the li terature, which is MRCl 138-262 
(| Nesvadba et al.l l2006b[ ) with densities in the range of 
240—570 cm~'^. The values in both galaxies are very similar, 
suggesting that densities of a few x 100 cm~'^ are likely typi- 
cal for the luminous optical emission line regions in HzRGs. 



Fig. 6. Relationship between electron temperature and 
[QUI] line ratios for an electron density of 500 cm~^. 
The green shaded area shows the temperature range for 
MRC0406-244 derived from a 3a hmit on [OIII]A4363, the 
red area shows the range of temperatures from the ratio 
based on the measured [OIII]A4363 flux in TXS0828-f 193, 
allowing for the 30% uncertainty. The inset shows a 
zoom on H7 and [OIIIJA4363 in the H-band spectrum 
of TXS0828-I-193. Blue lines indicate the expected emis- 
sion lines with the redshift and line width measured from 
[OIIIJA5007. The lower panel shows the fit rediduals. 



4.4. Electron Temperatures 

We detect the faint [OIIIJA4363 emission line in the H- 
band spectrum of TXS0828-I-193, which is the first time 
that this line has been measured in the spectrum of a 
z> 2 HzRG at > 3tT significance. We detected the line 
in a 5 X 5 pixel aperture in the H band centered on the 
bright southern knot. The spectrum around the wavelength 
of [OIIIJA4363 is shown in the inset of Fig. [Sj We use 
the ratio of emission line fluxes F of the 3 [QUI] lines 
at 4363A, 4959,A, and 5007A, respectively, F4363, F4959, 
and F5007, 15-5007+4959,4363 = (-^4959 + ^5007)7^4363 to esti- 
mate the electron temperature of TXS0828-H193. We find 
R.5007+4959,4363 ~ 173 ±52, Correcting for extinction. For an 
electron density Ug ~ 500 cm~'^, and allowing for ^ 30% 
uncertainty of the [OIIIJA4363 flux measurement, this cor- 
responds to an electron temperature of ^ 10100 — 12500 K 
(see Figure [6]). 

We do not detect [OIII]A4363 or H7 in MRC0406-244. 
This may be due to the somewhat lower redshift of this 



source, which causes both lines to fall at a wavelength 
where the transparency of the night sky and SINFONI's 
efflciency are not very favourable. Placing a 3ct limit 
on [OIII]A4363 in 1VIRC0406-244, we find an extinction- 
corrected R5007+4959.4363 = 55.6, which implies an electron 
temperature T°'^'^^ < 1.5 x 10"* K. 

Measuring electron temperatures from R[o///] is the 
most direct method, but obviously requires high-quality 
spectroscopy of high- redshift gala xies, and exposur e times 
that are often prohibitively large. iHumphrev et al.1 (|2008f ) 
explored alternative, temperature-dependent line ratios us- 
ing longslit spectra of rest-frame UV and optical emission 
lines. For 1VIRC0406-244 and TXS0828-f 193, which are both 
included in their sample, they find TO^oe ^ 11700^1^^!^ K 
and for several line ratios in TXS828-I-193, they find tem- 
peratures in the range Tf^^ = 14000 - 22000 K. While 
their results differ from ours if taken at face value, the dif- 
ferences are insignificant, given the large uncertainties of 
cither measurement. 
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4.5. Ionized gas masses 

We use the extinction corrected H a fluxes (ij4.2p to estimat e 
ionized gas masses. Similar to iNesvadba et al.l ()2006bD . 
we assume simple case B recombination, and use the 
Ue ~ 500 cm^'^ electro n density measured in MRC0406- 
244 (see also § 4.2 of INesvadba et al.l 1200 6b'. in partic- 
ular equation 1). Total Ha fluxes in MRC0406-244 and 
TXS0828-fl93 correspond to log [L/Lq]";!"^ = 45.2 erg 
s~^ and log [L/Lo]^^^^ = 45.3 erg s~^, respectively (or 
log[L/L0]°*°^ = 44.5 erg s'^ and log [L/Lgf^^^ = 44.4 
erg without extinction correction). Estimating the 
gas masses in each spatial resolution element, and sum- 
ming over the full emission of the two galaxies, we flnd 
log [AfH„/M0]""°^ ~ 10.6 and log[MHii /M^f^^^ - 10.7 
for MRC0406-244 and TXS0828-H193, respectively. All val- 
ues are also given in Table [71 

We caution that ionized gas mass estimates at z> 2 
can be precise to an order-of-magnitude level only. In par- 
ticular, extinction and electron densities involve measur- 
ing relatively faint lines, H/3 and [Sn]AA6716,6731, and the 
physical parameters are strongly non-linear function of the 
line ratios. In spite of signal-to- noise ratios ^ 10 in the faint 
lines, which is excellent compared to IFU studies of g alax- 
ies at high redshift (e.g. iForster Schreiber et aT]|2006f ) this 
will add uncertainties of factors > 2 in each estimate and 
for typical values measured in our sample. Correcting for 
the blend of the Ha+[NII] lines, we minimized the contri- 
bution of Ha in the flt, which will result in underestimating 
the Ha flux. Overall, we find that these uncertainties will 
amount in total to less than ±0.5 dex, safely placing our 
estimates into the range of a fewxl0^° Mq. 

A few studies have previously estimated ionized gas 
masses in HzRGs based on slit-spectra and narrow line 
images, alth o ugh w ith significantly larger uncertainties. 
lArmus et al.l (|l998l ) report on a giant nebulosity in 
4C-hl9.71 at z~ 3.6, which they found through [OIII]A5007 
narrow band imaging, estimating gas masses of ~ 10^~^ 
M0. From [OIIJA3727 longslit spectr oscopy of a sample o f 
52 HzRGs taken from the 3CR, iBaum fc McCarthy! (pOOl 
find that the ionized gas masses in HzRGs increase linearly 
with redshift, which may be related to larger quantities of 
cold gas in and around massive galaxies at high redshift. 
They also find that the 10 HzRGs at z > 1 have ionized gas 
masses of several IO^Mq, which is likely a lower limit, due 
to partial coverage of the source with the slit, and because 
[OIIJA3727 is susceptible to extinction. (In i )4.2l we found 
Ay = 1 — 4 mag for the extended emission line regions 
of our targets.) Comparing their result for MRC0406-244 
and ours, we can loosely estimate that correction factors 
of a few are required to compensate for systematic offsets. 
This suggests that the intrinsic masses of their sample will 
plausibly be of order lO^^M©. 



I Humphrey et al.l ()2008D measured Ha and H/3 fluxes for 
9 powerful radio galaxies during the "Quasar Era" with the 
ISAAC longslit spectrograph on the VLT. Unfortunately, 
H/3 is not or only marginally detected in most of their 
sample, making a direct comparison with the extinction- 
corrected gas masses of our sample difficult. However, for 
the uncorrected Ha fluxes, we find similar values within 
factors of a few as in our sample, which suggests ionized 
gas masses of order lO"'^" Mq may be common in HzRGs. 



5. The nature of the emission line region 

Rest-frame optical integral-field spectroscopy has been 
published for a few dozen gala xies at z^ 2 of dif- 
ferent types and m asses (e.g. iFo rster Schreiber et al.l 



200 1: INesvadba et al. "2007b". 2006a: Swinb ank et al 



La w et al.ll2007t INesvadba et al.,,2008; .Genzel et al 
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Wr ight et al.ll2007l: iBournaud et al.ll2008i ). including galax- 
ies undergoing mergers, candidates for relatively quiescent 
disk rotation and galaxies that are likely disks in forma- 
tion, and galaxies showing the signs of starburst-driven 
winds. The sample includes galaxies with and without 
AGN. Compared to these samples, our radio galaxies have 
several unique properties which allow us to constrain the 
mechanism which is driving the gas dynamics. We will 
briefly discuss the different possible scenarios. 

5.1. Merging 

If the two bubbles represented the interstellar medium of 
two massive galaxies in course of a merger, we should 
flnd a continuum peak in each bubble, representing the 
stellar component of each galaxy of the merging pair. 
However, we find only one continuum source in MRC0406- 
244 and TXS0828-I-193, respectively, which is in the center 
of each source, where the velocity of the ionized gas changes 
abruptly between the blueshifted and the redshifted bubble. 
(For MRC0316-257, we did not detect the continuum.) This 
does not rule out that these galaxies are related to a major 
merger, however, the two components must be closer than 
the resolution of our data, about 5 kpc, and the two bubbles 
cannot be gravitationally bound to one of the galaxies each. 
(Note that our continuum images are not deep enough to 
rule out minor mergrs.) 

Broad-band morphologies of HzRGs are in agree- 
ment with this interpretation. Broad-band morpholo- 
gies of HzRGs are known to be often irr egular (e.g., 
Ivan Breugel et al]|l998l : [P"entericci et al.ll200ll ). however, in 
many cases and depending on redshift, this may be due 
to emissi on line conta mination, as has been suspected be- 
fore (e.g. iRush et all 1 19971. and is now evident from the 
line and continuum morphologies extracted from our data 
cubes. A number of HzRGs, including galaxies with HST- 
NICMOS imaging, have very compact morphologies and el- 
liptical light profiles consistent with being merg er remnants 
(jvan Breugel et al.|[l998l : iPentericci et al.ll200l( ). 



5.2. Disk rotation 

Disk rotation and merging are difficult to distinguish at 
z~ 2 due to the relatively coarse spatial resolution of the 
data relative to the relevant size scales, typically several 
seconds of arc (see INesvadba et al.l l2008l for a more de- 
tailed discussion). However, all candidates for disk rotation 
at high redshift have more gradual velocity offsets. Unlike 
other high-redshift galaxy populations, the emission line 
regions of HzRGs are well resolved spatially, and should 
show a characteristic velocity pattern known as "spider di- 
agram" if their ki nematics were dominate d by disk rotation 
(see, e.g. Fig. 1 of lKronberger et al]|2007l ). However, the ob- 
served abrupt velocity offset within a spatially unresolved 
region characterized by strong line splitting (Fig. [7]) is at 
odds with galaxies showing disk rotation at low and at high 
redshift. In addition, the emission line morphologies are of- 
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Fig. 7. [OIII] emission lines extracted from a small aper- 
ture (0.375" X 0.375") near the center of MRC0406-244. Line 
splitting is evident in both lines, with three narrow compo- 
nents and an underlying, broad component. The mid com- 
ponent in [OIIIJA4959 is affected by a night sky line resid- 
ual. To guide the eye, we also show Gaussians for each com- 
ponent of the [OIII]A4959 and [OIIJA5007 line respectively. 
Note that these are not strictly fits, but have the same line 
width and redshift for both [OIII] lines. The red component 
of this line appears fainter than the corresponding compo- 
nent of the [OIIIJA5007 line, because of an underlying broad 
component of the [OIIIJA5007 line superimposed, which we 
ignored. 



ten irregular or filamentary, which would not be typical for 
rotating gaseous disks (see Fig. |8] to [10] ) . 

Moreover, the large velocity offsets and line widths sug- 
gest unplausibly large enclosed masses of often 10^'^"^'* 
M0 within few 10s of kpc (compared to s tellar masses 
of few xlO" Mq; ISevmour etld] |2007[ ). if the gas 
was to be bound gravitationally (see, e.g.. Table 1 of 
iBaum fc McCarthvl l20001. Fo r HzRGs tha t are s urrounded 
by galaxy overdensities, ,Venemans et al.l (|2007f ) estimate 
virial masses of few xlO^* Mq for the entire dark matter 
structure out to radii of 1 — 2 Mpc, based on observa- 
tions of the velocity dispersion of the companion galaxies. 
If the observed emission line kinematics of HzRGs were due 
to rotation, then a significant fraction of this mass would 
have to be within the inner 20-30 kpc of the structure, at 
odds with what we know abo ut dark matter halo mass pro- 
files. iNesvadba et al.l ()2006bD argued that the ionized gas of 
MRGl 138-262 has a kinetic energy which is about similar 
to the binding energy of a galaxy with a mass of few x 10"'^^ 
Mq, in which case much of the gas would not be gravi- 
tationally bound. We conclude that the observed velocity 
fields cannot be related to gravitationally driven motion, 
they are super-gravitational. 



5.3. Outflows or in fall? 

Directly measuring the direction of a gas flow in extra- 
galactic objects is inherently difficult, in particular at high 
redshift. Observed velocity gradients are near the escape 
velo city of the dark halo s urrounding a massive galaxy 
(see INesvadba et al.l l2006bl for a more detailed discus- 
sion), which suggests that we may be witnessing the ion- 
ized gas being unbound from the host galaxy. However, 



in the next section we will argue that the gas dynamics 
are most likely related to the expanding cocoon of hot, 
overpressurized gas surrounding the radio jet in agreement 
with previous stti d ies based on longslit spect r oscopy (e.g., 
Clark et all 119981: IVillar-Martm et al.l 119991 : llnskip et al] 



20021: iBest et al.lll997D . 



Such an interaction will result in complex kinematics, 
in agreement with the observed large line widths. The co- 
coon likely forms from backflowing gas which has previ- 
ously been he ated by shocks at th e work ing surface of the 
jet. Although [Alexander fc PoolevI (|l996( l pointed out that 
this gas will likely not move far from the hotspot and may 
mostly contribute to inflate the cocoon, this makes a more 
detailed assessment of the flow direction of the ionized gas 
worthwhile. 

One way to constrain this direction is through a care- 
ful study of the radio spe ctral properties of th e two 
lobes in each source (see a lso iHumphrev et aI1l2008l ). (e.g., 
[Alexander fc Leahvl 119871 ). Radio emission from a source 
embedded in a magnetized plasma will be depolarized by 
a degree that depends on the column of the transversed 
medium. Thus, the near side of the radio source will ap- 
pear more strongly polarized than the far side (the " Laing - 
Garrington effect"; [LaingI [19881 : [Garrington et all [l988( i. 
Moreover, the observed brightness of the radio hotspots 
on the near side of the ra dio source may be enhance d due 
to Doppler-boosting (e.g.. [Alexander fc Leahyiri987D . The 
hotspots have younger ages and accordingly flatter spectral 
indices than the surrounding lobes. Thus, at high frequen- 
cies, the relative contribution of the hotspot will contribute 
more strongly to the overall emission from the lobe, and 
the spectral index on the near side of the radio source will 
appear flatter than that on the far side^ 

Comparing with the observations of lCarilli et al.l ([1997[ ). 
we flnd in each of our galaxies that the northern radio lobe 
is more strongly polarized at 4.7 GHz and at 8.2 GHz than 
the southern lobe. Moreover, the northern lobe has a some- 
what flatter spectral index in all galaxies. In all three cases, 
this is also the side of the blueshifted bubble of ionized gas. 
This strongly advocates that the blueshift represents a net 
outflow relative to the central galaxy. 

Outflows of gaseous material are also ob served in power- 
ful rad io galaxies at lower redshifts, z< O.l. lMorganti et al] 
([2005[ ) find broad, blueshifted HI absorption lines at rela- 
tive velocities of > 1000 km s^^ from the systemic redshift 
in 7 out of 8 sources. This is consistent with AGN-driven 
outflows of neutral gas, corresponding to mass loss rates of 
up to ~ 5OAf0 yr"^ 

With the simple premise that the extended emission line 
regions in HzRGs represent energetic outflows, we will now 
constrain the basic dynamics of the gas, before arguing in 
the next section, that the outflows are most likely related 
to the radio jet. 

5.4. Dynamical Time Scales 

Emission line regions extend over fairly similar radii in 
all three galaxies, 10 kpc x9 kpc, 15 kpc x9 kpc, and 
14 kpc x9 kpc for MRC0316-257, MRC0406-244, and 
TXS0828+193, respectively (see Table [6]). We can use the 
size of the outflows along the axis of the radio jet and the 
relative velocity between the blueshifted and the redshifted 
bubble in each galaxy to estimate the crossing time. Using 
the velocities given in Table [HI we find dynamical timescales 
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of Tion ~ few X 10^ yrs, which is similar to ages inferred for 
radio jets (e.g., iKaiser fc AlexandeJfT99l . See Tabled for 
the ages derived for individual galaxies. These estimates 
can only be accurate within factors of a few, given the un- 
known radial density profile of the gas prior to the radio- 
loud phase, the effects of projection and relatively coarse 
spatial resolution of our data. 

Alternatively, we can assume that the velocity gradi- 
ents observed in individual bubbles are dominated by the 
lateral expansion of the bubbles, as observed in low-redshift 
narrow-line regions being inflated by interactions with the 
overpressurized, exp anding cocoon of the radio jet (e.g., 
ICapetti et al.l[l999( ). HzRGs h a.ve very irregular environ - 



ments with gaseous halos (e.g., IVillar-Martm et al 



satellite galaxy overdensities (e.g. IVenemans et al 



2003 ) 



20071) 



and perhaps infalling gas clouds, which will influence the 
distribution and density of the sourrounding medium. This 
will make it difficult to reliably measure the relatively small 
velocity gradients within an individual bubble. However, if 
we only use clouds with relatively smooth internal veloc- 
ity gradients, in particular the blue (northern) bubbles of 
MRC0316-247 or of MRC0406-244, we find internal velocity 
gradients of ~ 120 — 300 km s~^, corresponding to dynam- 
ical time scales of a few xlO^ yrs, which is consistent with 
the dynamical time scales estimated from the velocity gra- 
dients along the radio jet axis. 

5.5. Kinetic Energies 

In §5.2 and §5.3 of iNesvadba et al.l ()2006bD . we applied two 
methods to estimate the injection of kinetic energy into the 
outflows: 

(1) A strict, but presumably not very precise lower limit, 
based on the observed ionized gas mass, the bulk motion, 
V, and the turbulent motion, a, i.e., E = l/2m(w^ + cr^). 

Velocities and line widths are taken from the kinematic 
maps shown in Figures [51 to [TUl and thus represent the local 
values for each spatial resolution element. We discarded 
the very broad lines near the center of the galaxy, which we 
believe are largely due to overlaps between the red and the 
blue bubble. Within each bubble, line widths are larger than 
the velocity changes, so that they give a fair representation 
of the intrinsic line widths within the spatial resolution of 
our data. 

Since we neglected the likely presence of other gas 
phases, as well as uncertainties related to projection, this 
can only be a conservative lower limit. We do not assume 
pressure equilibrium, but that velocities are dominated by 
ram pressure. 

(2) An energy estimate based on the assumption that the 
ionized gas is being entrained and accelerated through adi- 
abatic inflation of a bub ble of hot o verpressurized ma- 
terial, which follows, e.g., iDvson fc Williams (1980), and 
which yields E = 1.5 x lO^^r^gU^^QQQno.s erg s~^. rio indi- 
cates the size of the emission line nebulae in units of 10 
kpc, uiooo the velocity in units of 1000 km s~^, and no.5 
the ambient gas density of the warm interstellar medium 
prior to being swept up by the bubble in multiples of 0.5 
cm^"^. This approach is highly idealized given the com- 
plexity of the turbulence and bulk motion observed and 
expected from hydrodyn amical simulations of radio jets 
(e.g.,[Kr ausc 2005; Saxto net al.ll2005l : iGaibler et aLll2007t 
ISutheriand fc BicknelL2007i r but may serve as a rough an- 



alytical surrogate to more detailed estimates. We used the 
velocity maps shown in Figures [5] to [TO] to estimate the to- 
tal velocity offsets between bubbles, and set the expansion 
velocity in each bubble to half the total value. Values for 
individual galaxies are given in Table [51 

Although our estimates of velocities and line widths are 
precise to 10%-20%, each approach are based on a set of 
very simple assumptions. While this is appropriate given 
the observational difficulties of high-redshift studies, it does 
introduce astrophysical uncertainties that will make each 
estimate precise to factors of a few only. 

Applying method (1) to our present data set, we find en- 
ergy injection rates of 3 x 10^^ erg s~^ and 1.3 x 10^** erg s^^ 
for TXS0828-hl93 and MRC0406-244, respectively, for mass 
estimates without extinction correction. This represents a 
very conservative lower limit on the energy injection rate. If 
using extinction-corrected mass estimates instead, we find 
energy injection rates of ~ 2 x 10^^ erg s~^ are required for 
both galaxies with ionized mass estimates, MRC0406-244 
and TXS0828-f 193. (Individual values are given in Tabled) 

For MRC0316-257, we do not have an explicit estimate 
of the ionized gas mass, because at the redshift of the 
source, z~ 3.1, Ha cannot be observed from the ground, 
and H/3 is very faint. The velocity maps shown in Figures [H 
to [TOl , however, suggest overall very similar properties. 
Under the premise, that the ionized gas mass will not 
greatly differ, we would hence expect to also find simi- 
lar estimates of the kinetic energy injection rate. For the 
10 galaxies in the sample of iBaum" fc McCarthvl (l2000D at 
z> 1, we find very similar results. Using their measure- 
ments of the ionized gas masses, velocities, and line widths 
and our method 1, we find total energy injections of 10^'^"^^ 
erg over dynamical timescales of about 10*" yrs. Note that 
this is a strict lower limit. If we correct for extinction as 
found from our SINFONI data set, and uncertainties in the 
alignment between the longslit and the direction of largest 
velocity shear, these estimates will easily increase by about 
an order of magnitude. This certainly highlights the ubiq- 
uity of powerful outffows of ionized gas in HzRGs. 

Method (2) does not depend on the ionized gas mass, 
and can therefore be applied to all three galaxies. We used 
the velocity maps shown in Fig. 8 to estimate the total 
velocity offsets between bubbles, and set the expansion ve- 
locity in each bubble to half the total value. Values for 
individual galaxies are given in Table 6. 

We find that energy injection rates of i? ^ 0.7 — 4 x 10^^ 
erg s^^ are required to power the outflows (see Table [HI for 
individual values). Note that velocities are not corrected for 
orientation effects. Assuming that inclinations will typically 
be around 45° or higher, intrinsic velocities may be factors 
~ 2 — 3 higher, which would result in kinetic energies of up 
to several xlO**^ erg s^^ and total kinetic energy injections 
during an AGN lifetime of order 10^° erg. 



6. The role of the radio jet 

6.1. Kinematics 

Our study provides direct evidence for outflows of ionized 
gas driven by p owerful AGN during the "Quasar Era". 
INesvadba et all (|2006b r) argued in the case of MRC1138- 
262 that the energy injection, mass outflow rates, and 
dynamical time scales satisfy the criteria set by observa- 
tions of massive galaxies at low redshift. In particular, the 
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outflow time scales are short enough to explain the ob- 
served relative enhancement of a elements over iron as a 
function of velocity dispersion in early-type galaxies (e.g., 
iPipino k Matteucci|[20M ; iPipino et al.ll2007D . Direct stud- 
ies of galaxies undergoing phases of strong AGN feedback 
like HzRGs are ideal to investigate by which mechanism the 
energy output of the AGN is being converted into thermal 
and kinetic energy of the surrounding cold gas. Broadly 
speaking, two groups of processes have been proposed in 
the literature (1) Outflows driv en by the radiative e nergy 
output of the AGN (e.g.. iKind 12001 ISazonov et al.l l2005t 
iMurrav et al] 120051 : ICiotti fc Ostrikeil I2007D . (2) Outflows 
related to mechanical interactions with the relativistic, 
synchrotron-emitting plasma of radio-loud AGN. These 
were mai nly proposed to ove rcome the cooling- flow problem 
(e.g .[Brn ncy & Tabor 19951: iRuszkowski fc Begelmanll2002t 
iKaiser fc Binnevll2OO30 . but are also addres sing the impact 
of th e jets in indiv idual radio galaxies (e.g., Mellema et ahl 
20021: iFragile et al . 2004; Sa xton et al.ll2005l: iKrausd 120051: 



Sutherland fc Bicknel]||2007l: IKrause fc Alexanderj|2007l)" 



While both classes rely on physically plausible mecha- 
nisms, their actual impact on the surrounding gas depends 
strongly on the efficiency by which the energy output of the 
AGN is being converted into thermal and kinetic energy of 
the surrounding gas. This is best addressed through direct 
observations of systems that are in a phase dominated by 
powerful, AGN driven outflows, like the galaxies of our sam- 
ple. Overall our analysis favors a picture where the radio 
jet plays a dominant role in driving the gas, based on ar- 
guments related to the observed geometry, timescales, and 
kinetic energy of the flows: 

Orientation relative to the jet axis: In the 3 galaxies 
of our sample, the extended emission line regions are elon- 
gated along the axis of the radio jets, and are well aligned 
with the jet axis. The same is suggested by the overall ve- 
locity field of the galaxies - we clearly distinguish a red- 
shifted and a blueshifted bubble, suggesting a back-to-back 
flow, with a sharp transition near the center of the galaxy. 
The size of the transition regions appears dominated by the 
spatial resolution of our seeing- limited data. In the previous 
section, we also showed that the blueshifted gas appears to 
be aligned with the near side of the radio source in all three 
cases. This is a prerequisite, if the ionized gas is being ac- 
celerated away from the galaxy by the expanding cocoon of 
hot gas surrounding the radio jet. 

Spatial Extent: The region of high surface brightness 
emission extends to radii that are smaller than the size of 
the rad io lobes. Extending this argument, iNesvadba et al.l 
(|2007af) find that galaxies with compact radio sources do not 
have spatially extended line emission, but kinetic properties 
suggesting that the gas kinematics is driven by the same 
underlying proce sses. In agreement with spectral studies 
of ra dio jets (e.g. iMurgia et al.l 1X9991 : lOwsianik fc Conwavl 
Il998i) . this is a natural result if compact radio sources 
are predominantly young. For similar entrainment rates 
as found in our sample, and typical ages of a fewxlO^"^ 
yrs estimated for compact radio sources, these galaxies 
should not yet have entrained substantial fractions of the 
interstellar medium, i n agreement with observations (see 
INesvadba et al.ll2007al f or details). 



Villar-Martm et all ((200l . [2003L [200l [200l find faint 



rest-frame UV emission surrounding HzRGs, with line 
widths and velocity gradients that suggest that the gas 
does not participate in the violent processes in the inner 



halo, that we are tracing with our rest-frame optical data. 
These kinematically quiescent outer halos are always found 
at radii beyond the radio jets, which is highly suggestive in 
a scenario where the outflows are dominated by interactions 
with the radio jet. 

Time scales: In ij5.41 we estimated dynamical time scales 
of a few xlO^ yrs from the characteristic velocities and 
spatial extent of the flows. In spite of large uncertainties, 
this is in the typical range of age s estimated for radio 
jets fe.g.. 'Kaiser fc Alex ander Il997t [Kaiser fc Best! 120071: 
Blundell & Rawlings 199^. 

Energies: Anticipating the results of %M we estimate that 
about 10% of the kinetic energy flux of the jet is trans- 
formed into kinetic energy of the gas. This is fully consis- 
tent with a picture where the radio jet is the primary driver 
of the observed gas kinematics. 

These large coupling efficiencies may be somewhat sur- 
prising. Early models of radio jets indicated that powerful 
jets like the ones in our sample (in the energy range of FR 
II radio galaxies), will pass through the ambient medium 
rather undisturbed, while low-power jets (in particular FR I 
galaxies) may deposit most of their ener gy in the surroun d- 
ing medium after being disrupted fe.g.. lDe Yoim3ll99lD . 

More recent models with better resolution, and larger 
density contrast between jet plasma, ambient medium, and 
embedded clouds, suggest a more subtle picture. Large 
jet kinetic energies of few xlO'*^ erg s~^ and highly in- 
homogeneous ambient media are clearly required to ad- 
equately describe powerful radio galaxies at z~ 2, and 
several recent models have now investigated t h is (e.g. 
'Saxton et al.' 2005; 'Sutherla nd fc Bicknelll l2007t iKrausel 
2005; Krause & Alexander 2007^. The basic properties of 
the simulated cocoons are not unrealistic compared to what 
is observed, e.g., internal velocities of order 1000 km or 
aspect ratios o f ~ 3 between the radial and the lateral axis 
of the cocoon (|Saxton et all [20051 : lKrausd[2005h . 

However, the model parameters are adapted to lower 
redshift galaxies with somewhat lower jet eneriges, and 
lower fractions of cold gas. In particular, the environ- 
ments of HzRGs suggest the galaxies are surrounded 
by halos of cold gas with complex substructure, rather 

than well-settled dis ks of cold gas as adopted by 

[Sutherland fc Bicknelll ([200l . 

ISaxton et alj (|2005|) point out that a clumpy environ- 
ment may destabilize or even disrupt the jet by scattering 
on several clouds, and through Kelvin-Helmholtz instabili- 
ties induced by the overall level of turbulence rather than 
individual jet-cloud interactions. It is plausible that this 
will enhance the efficiency with which the kinetic jet en- 
ergy is being deposited into the cold gas even for the most 
powerful jets. Models better adapted to the properties of 
powerful radio galaxies at Z'^ 2 may be an interesting way 
of investigating this further, including the kinematic con- 
straints now available from integral-field spectroscopy. 



6.2. Multiphase gas content of HzRGs 

6.3. Mass budgets 

The extended emission line regions of HzRGs do not only 
have extreme kinematic properties, but also an unusually 
large content of ioriized g as, few xlO^*^ Mq ( il4.5i see also 
iBaum fc McCarthvll2000r ) . This exceeds the amount of ion- 



12 



Nesvadba et al.: Jet-driven AGN winds in Radio Galaxies in the "Quasar Era" 



ized gas in any other high-redshift galaxy population dis- 
cussed in the literature by several orders of magnitude. 

For example, the same mass estimate applied to the 
total Ha emission of 2 act i vely star-forming "BX"- 
galaxies measured by lErb et alj (|2006f ). would only corre- 
spond to a few xlO^ Mq, and most of this gas is proba- 
bly related to s tar-forming regions and does not participate 
in an outflow. iNesvadba etahl (|2007bD used SINFONI to 
obtain spatially-resolved data of the starburst-driven wind 
in a submillimeter selected, strongly star-forming galaxy 
(SMG), SMMJ1401H-0252. They find that in spite of high 
star- formation rates of SFR^^ 300Mq yr~^, only ~ 10% of 
the total Ha flux is in a blueshifted component - the wind 
- which would correspond to ionized gas masses of a few 

XlO^ Mq. 

Overall, the stellar masses and star-formation rates in 
z~ 2 HzRGs closely resemble those of SMGs at similar 
redshifts. Moreover, several IO^^Mq of cold molecular gas 
traced by CO millimeter line emission have been observed 
in SMGs (e.g., jNerTeraL 2003; Grcvc ct al.'200E |) as well as 



in HzRGs (e.g.TlPapadopoulos et al.. 2000 : De Breuck et al.l 
[20031 [200a iKlamer et all 120051) . Large masses of cold gas 



are a necessary prerequisite to fuel a strong star burst. 
However, not all HzRGs have been detected in CO. For 
TX0828-t-193 specifically, which appears to have several 
10^° M© in ionized gas(> j4.5p . Nesvadba et al. (in prep.) 
did not detect CO line emission, but place an upper limit 
of ^ 1 X 10^°Mq in cold molecular gas within the galaxy 
(for a ULIRG-like conversion factor from CO luminosity 
to H2 mass, similar to that adopted for other high-redshift 
CO-luminous galaxies). Taken at face value, this is less than 
the amount of ionized gas we observe with SINFONI. Thus, 
it appears that the ratios of ionized to molecular gas in 
HzRGs, are of order unity or greater, Rf^^^i > 0{1), com- 



pared to i?- 



SFG 
ion.mol 



< )(io 



-(3-4) 



) for actively star-forming 



galaxies at z 2 (SFGs) including SMGs. (Either value is 
affected by observational uncertainties of factors of a few, 
but the large contrast between star-forming galaxies and 
radio galaxies makes this comparison robust nonetheless.) 
This certainly highlights the importance of the observed 
outflows for the evolution of the host galaxy, and it is in- 
teresting in this context that some HzRGs show the signs of 
strong star-formation (like large reservoirs of molecular gas, 
strong PAH and submillimet er emissio n, 'Archibal d et al.l 
[2001; Rculand ct al. 2004; Sc vmour et a l. 2008), while oth- 
ers do not. This is certainly to be expected for a population 
that marks the transition from an actively star-forming to 
a passively evolving evolutionary stage. A detailed study of 
the multiphase gas content of HzRGs for a more represen- 
tative sample will be crucial to assess this question directly. 



6.4. Ionizing mechanism 

Closely related is the question which mechanism domi- 
nates the ionization of the gas? This is particularly im- 
portant since our gas mass estimates (similar to those 
of iBaum fc McCarthy! I2OOO0 do assume that most of 
the gas is photoionized, and in fact, the bolometric en- 
ergy output of HzRGs seems sufficient t o postulate that 
the n ebulae may be matter bounded ([Nesvadba et al.l 
l2006bf) . Photoionizat ion is also favoured by the mea- 
sured UV line ratios ([Villar-Martm et al.|[T999f) . However, 
the gas kinematics suggests that shocks play a domi- 



nant role for the overall energy budget of the outflows, 
and hence it would appear natural that shocks would 
be an important ionizing source (a well-known puzzle , 
for more detailed discussions see e.g. iClark et al.l [l99l 



IVillar-Martin et all [19971 [Allen et al 



et al.l 
1995L 



200s ). Using 
1996li for the 



the 

models of iDopita &: Sutherland! (|1995l I1996D for the line 

emission of fast shocks, we find that by assuming that shock 
ionization dominates, gas mass estimates would be about 
an order of magnitude larg er than what we estimat ed from 
case B recombination (see [Nesvadba et al.l [2006bl for de- 
tails). In this sense, and in spite of uncertainties of factors 
of a few, our gas mass estimates are rather conservative. 

Ultimately, photoionization may explain why we ob- 
serve larger amounts of ionized gas by several orders of 
magnitude than what is found in low-redshift radio galax- 
ies, where fast, AGN-triggered outflows are predominantly 
seen in neutral gas through broad, blueshifted (^ 1000 
km s"^) HI absorption fines ([Morganti et all [20051) . This 
dichotomy may be mainly due to the different radiative 
properties of the AGN rather than fundamentally differ- 
ent acceleration mechanisms for the outflows; HzRGs are 
harboring luminous, obscured quasars, and are embedded 
in large quantities of cold gas, whereas low-redshift radio 
galaxies seem to be radiatively inefficient and reside in hot. 
X-ray emitting halos. If this is the case, then we may sig- 
nificantly underestimate the frequency of AGN driven out- 
flow s, simply because the y are observationally more subtle 
f see [Morganti et an[2005l for a more detailed discussion). 



6.5. Entrainment rates 

The amount of ionized gas and the dynamical time scales 
provide a rough estimate of the entrainment rate, as the 
cocoon plasma interacts with the ambient medium. The ef- 
ficiency of the overpressurized plasma in entraining colder 
(T< 10^ K) material is one of the crucial parameters in 
characterizing the efficiency by which the radio jet inter- 
acts with the ISM of a galaxy, and is a parameter that is 
now being addressed by detailed hydrodynamic simulations 
fe.g.. [Krausj[2005t [Sutherland &: Bicknel][[2007[ ). However, 
since the amount of entrained material will strongly depend 
on the "macro" and "micro" -physics of the jet and ambient 
material, including several orders of magnitude in pressure, 
density, and size scales, it is difficult to derive this quantity 
from first principles. 

The ionized gas we see will most likely be on the surfac e 
of filaments with low filling factors ([Nesvadba et al.[[2006b[ ). 
that are being swept up, accelerated, and ionized by the hot 
medium surrounding the radio jet. We can therefore give a 
purely empirical estimate of the entrainment rate, assuming 
that it roughly corresponds to the outflow rate of ionized 
gas, Mion, and the dynamical timescale, r. 

Recasting the discussion in [Nesvadba et al.) ([2007a[ ). we 
approximate the rate of entrained material, M(,ntr,9 (in mul- 
tiples of IO^Mq), from the ionized gas masses of xlO^° Mq 
and dynamical timescales of few xlO^ yrs. This yields a 
simple scahng, Mcntr.g = fc fp,fgcomMion tj, where T7 indi- 
cates the dynamical time scale in units of 10^ yrs. fc, fp^, 
and fgeorn parametrize our ignorance of the jet advance 
speed, radial density profile, and geometry, respectively. If 
we assume that all are of order unity for the relevant radii 
and time scales, we find an average entrainment rate of 
Mentr ~ lO^~^M0yr~^. The cylindrical shape of the out- 
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flows suggests that the outflows do cover large solid an- 
gles, since they are related mainly to the expansion of the 
cocoon, and contrary to what may be expected from the 
pencil-like morphology of the jet. 

The large ratio of ionized to molecular gas of order unity 
will make this estimate relatively robust. Nonetheless it will 
be a lower limit, given that we have no direct constraints 
on the amount of hot gas in the cocoon, and CO only 
traces the col d molecular ga s with temperatures of order 
~ 10 - 100 K. iKrausd (|2005[ ) argue that the mass of X-ray 
emitting gas within the cocoon may not great ly exceed the 
mass e stimates for the colder gas given above. iCariUi et al.l 
(|2002D estimate that the X-ray emission of MRCl 138-262 
may originate from up to ^ 10^^ Mq of hot gas, however, 
more recent X-ray studies of HzRGs favor the view that 
much of the emission originates from Compton scattering 
of photons from the cosmic microwa ve background on th e 
relativistic electrons of the jet (e.g. lOverzier" et al.ll2005D . 
and accrodingly less X-ray emitting gas. 

6. 6. Coupling Efficiencies 

If jets are to cause the observed outflows, then they must 
provide at least the kinetic energy observed in the gas. The 
synchrotron emissivity of a radio source is a strong function 
of the ambient medium into which a jet expands. It is there- 
fore difficult to estimate the intrinsic kinetic energy of the 
relativistic plasma from the luminosity of the radio source. 
Several methods have been prop osed in t he literature (e.g . 
Bicknell et a l* ^199?; Willott ct al. 1991; IWan et all [200a 
Birzan et al.li2004; Kaiser fc Best 2007.) which suggest that 
of order 0.1%-10% o f the kinetic energy of a j et is emitted in 
the radio. Similar to lNesvadba et al] (|2006bD we adopt two 
methods to estimate the jet kinetic energies of the galaxies 
in our sample. (1) A simple scaling, assuming that 10% of 
the kinetic energy is radiated at frequencies between 0.1-1 
G Hz in the rest-fram e. (2) The empirical relationship found 
bv lWan et all ()2000[) between the observed flux at 178 MHz 
frequency and the kinetic radio power. With the two meth- 
ods, we find that the jets eject several 10^^ erg s~^ in kinetic 
energy. Given that the adopted values are close to the lower 
bound of conversion factors from observed to intrinsic lu- 
minosities, these values serve effectively as lower bounds. 
For the estimates of individual galaxies see Table El 

A critical quantity to gauge the overall impact of jet- 
driven AGN outfiows is the efficiency with which the jet 
kinetic energy, Ekinjet, is being translated into kinetic en- 
ergy of the outflow, i?fei„_gas, the coupling efficiency rj = 
Ekin,gas/Ekin.jet- Obviously, such an estimate will not be 
very precise, given the large uncertainties in each measure- 
ment. Nonetheless, we feel that quantifying this value at 
order-of-magnitude precision will already be useful. With 
the values given in Table [6l we find overall coupling effi- 
cie ncies of ^ 5 — 80%. 

I Kaiser fc Alexander! (|1997f) pointed out that for approx- 
imately adiabatic expansion, a radio source will become 
fainter as it expands. However, in our simple attempt to 
estimate jet kinetic energies, we did not take this into ac- 
count. This may explain why we find the highest efficiencies 
in TXS0828H-193, which is the largest source in our sample, 
by a factor ~ 2. We may simply underestimate the kinetic 
luminosity of the radio jet by factors of a few. We therefore 
favor a coupling efficiency of order 10%. Similar efficiencies 
have been found in radio galaxies with outflows of neutral 



gas at significantly lower redshifts (jMorganti et al.l2005f ). as 
well as in the X-ray cavities of massive low-reds hift galaxy 
clusters (e.g-TMcNama ra & Nulscn 2007; Birzan et al.ll2004l : 
'Ra ffertv et al.ll2006D . This may be a first indication that the 
basic mechanism in all these environments may not be fun- 
damentally different although they act on different scales. 

7. The global role of outflows in HzRGs 

One of the main questions we address in this study is on 
the frequency of highly turbulent emission line gas and en- 
ergetic outflows in HzRGs, similar to what we observed in 
our SINFONI study of MRC1138-262. Such outflows may 
be an important mechanism in terminating the vigorous 
starbursts observed in galaxies near the upper end of the 
galaxy mass function during the "Quasar Era" , if they are 
powerful enough to efficiently heat and unbind a signifi- 
cant fraction of the interstellar medium of a massive, gas- 
rich galaxy. The presence of such outflows in another three 
HzRGs with extended radio sources (among three galaxies 
observed) suggest that AGN-driven outflows are a common 
phenomenon in z~ 2 — 3 HzRGs. 

Overall, the physical conditions and gas kinematics in 
MRC1138-262 and the three galaxies of the present sam- 
ple are very similar. All four galaxies have typical line 
widths of FWHM - 800 - 1000 km s"^ over large areas, 
suggesting turbulent energy injection rates of ~ 10'''' erg 
s~'. Areas with larger widths may be due to overlap of 
adjacent bubbles. Ha fluxes indicate fairly similar ionized 
gas masses in all four galaxies. The physical conditions 
within the outflows seem also largely similar, with elec- 
tron densities of several hundred cm~'^ for MRC0406-244 
and MRC1138-262. Both TXS0828-hl93 and MRC0406- 
244 have [OIII] emission line ratios consistent with tem- 
peratures of ~ lO'' K. Relative velocities between indi- 
vidual bubbles are ~ 1000 km s~' in all cases, suggest- 
ing velocities near the escape velocity of a massive dark 
matter halo. Differences are mainly in the emission line 
morphology where MRC1138-262 seems particularly irreg- 
ular. The reason may be different evolutionary stages or 
perhaps that MRC1138-262 has a more extreme environ- 
ment. This is suggested by the particul arly large veloc- 
ity dispersion of companion galaxies (e.g.. ^^nemans et al.l 
2007; Milcy et al. 2006) and an exceptionally large rota- 
tion measure (jPentericci et al.ll2000l ). However, we cannot 
rule out that the greater complexity of MRC1138-262 is 
simply the result of different orientations, with the AGN 
of MRCl 138-262 being less inclined relative to the line of 
sight. MRCl 138-262 hosts an ob scured quasar evidenced by 
broad nuclear Ha line emission (|Nesvadba et al.ll2006b[ ). 

The global dynamical properties of the emission line 
gas however do appear very similar in all 4 galaxies of our 
current sample, in spite of this range in the appearance 
of individual sources. This suggests that powerful outflows 
in HzRGs likely are a common phenomenon. Studying the 
outflows in a larger sample of HzRGs with and without 
broad Ha lines may help resolving this issue, and we have 
now started observations of a significantly larger sample of 
~ 30 HzRGs with SINFONI. 

7.1. Subsequent Mass Assembly of HzRGs 

The properties of HzRGs suggest that they are near the 
end of their main epoch of active star-formation, as seen. 
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e.g., from the evolution towards very regular, probably 
elliptical con tinuum profiles in HzRGs at z~ 2 (e.g., 
Ivan Breugel et al. f998). This is certainly consistent with 
postulating that the AGN winds we observe are truncat- 
ing the star-formation epoch in the host galaxy. However, 
will the fingerprints of this epoch be preserved until z—0? 
HzRGs reside in overdense environments, and may dou- 
ble their stellar mass through subsequent sa tellite accretion 
(as p redicted by galaxy evolution models, iDe Lucia et all 

mm- 

We can roughly estimate the accretion-dominated 
growth of HzRGs from the number of compan- 
ion galaxies observed around many HzI lGs, including 
MRC 316-257 and MRC113 8- 262 (e s. iLe Fevre et al.l 



_ _ .Venemans et all l2007l l2005t ' iKurk et all l2004t 
Best et aTT 2003). Continuum luminosities and Ha line 



19961 



widths (e.g.. lKurk eralll2004[ ) imply masses of typically few 
xlO^M0. With up to ~ 50 companions, the central galaxy 
may gain of order 10^^ Mq, if all companions are accreted, 
correspoding about half its stellar mass. For gas fractions 
of order 50%, satellite accretion may fuel the formation of 
~ 20 — 30% of the total stellar mass in the descendent of a 
HzRG at z~ 0. 

However, this will not inevitably dilute the fingerprints 
of previous AGN feedback, such as the low content of 
cold and warm gas or truncated star-f ormation histories of 
mass ive galaxies at low redshift (e.g., iPipino fc Matteuccil 
l2004f ). Dynamical interactions in the dense environment of 
HzRGs will lead to premature ageing of the stellar pop- 
ulation in the satellites and effi cient gas removal , similar 
to low-redshift galaxy clusters. .Boselli et al.l (|2008l ) find for 
dwarf galaxies in the Virgo cluster, that ram pressure strip- 
ping alone may turn strongly star-forming irregular galaxies 
into gas-poor ellipticals within 1—2 crossing times (few 100 
Myrs for a HzRG halo) and lift the satellite onto the red 
sequence. (Interestingly, the red sequence se ems to fall in 
place at z~ 2, as shown [Kodama et al.ll200"7l for the over- 
densities surrounding MRC1138-262 and MRC0316-257.) 
The metal enrichment of the ISM will be particularly ef- 
ficient, resulting in a metallicity increas e of about 0.5 dex 
relative to the field (jBoselli et al.ll2008[). Comparing with 
the mass-metallicity relationship of iTremonti et all (120041 ) 
we find that this will compensate for the metallicity differ- 
ence between a 5 x 10^ Mq galaxy and a 10^^ M0 galaxy. 

AGN outflows may enhance these processes, either 
throu gh gas stripping when a satellit e crosses the out- 
flow llrwin et all Il987t iFuiital I2008D . or by increas- 
ing the e ntropy and pressure in the surroundi ng hot 
medium (|Nath fc Rovchowdhurvl 120021 : [McCart hv et all 
|2008[) . where an AGN-driven pre-heating phase at z~ 1 — 3 
appears sufficient to explain the observed temperature pro- 
files of massive galaxy clusters. In either case, the satellite 
galaxies will likely undergo "dry" mergers with low dissi- 
pation, and may have a surprisingly small impact on the 
properties of the low-redshift descendents of HzRGs. Dry 
merging may also explain why HzRGs seem to have smaller 
effective radii than massive spheroidal galaxies at low red- 
shift, as well as th e extended stellar envelopes of cD galaxies 
(|Schombertlll987[ ). 



7.2. The cosmological significance of jet-driven outflows 

With only 6 HzRGs with rest-frame optical IFU data, our 
sample is relatively small. However, the great similarity of 



the ionized gas masses and kinetic properties is intrigu- 
ing, and raises confidence that similar outfiows are com- 
mon among powerful radio galaxies during the era of lu- 
minous quasars. Moreover, large ionized gas masses, line 
widths, and velocity gradients have been reported previ- 
ously, based on observations with longslit spectrographs 
of radio galaxies at rcdshifts z~l and in a few cases 
below (e.g Baum & McC arthy! I2OOOI: IVillar- Mart m et all 
20021. 120031 il999: McCar tlwTt al.l[l996t iTa dhunter 11991 
Inskip et al.lr2002 : .Clark et al.lll998l : iBest ~al..,1997l l. 
In t j6.6l we argued that of order 10% of the jet kinetic en- 



ergy is transformed into kinetic energy of the gas. We will in 
the following assume that this is a typical value for powerful 
radio sources at z~ 2. We can then use the observed radio 
luminosity function of steep-spectrum radio sources to esti- 
mate the global energy release of jet-driven AGN feedback 
into the intergalactic medium at high redshift^ 

W e use the radio luminosity function of IWillott et al.l 
(|200lD . who find that the radio-loud population can be split 
into two classes, a low-luminosity population associated 
with galaxies with low-luminosity emission lines includ- 
ing FRI and FRII radio structures, and a high-luminosity 
population consisting of radio galaxies and quasars with 
strong emission lines, and with nearly exclusively FRII ra- 
dio sources. We only include this luminous population into 
our estimates, and conservatively use their set of model 
parameters which minimizes our estimate. The comoving 
number density of this population rises by nearly 3 orders 
of magnitude between z=0 and z=2, and has properties 
that generally agree very well with the sources we observe. 

We computed the luminosity function for redshifts z— 
— 4 in bins of z= 0.1, and use the peak of the co- moving 
number density at a radio-power at 151 MHz of £151 = 
10^^ W Hz~^ sr~^, which is sufficiently close to the ra- 
dio power of the sources in our sample. We then compute 
the co-moving number density per redshift bin, ^modi con- 
verting to our cosmology. Since the radio-loud phase, Tj^t 
is significantly shorter than the cosmic time within these 
bins, we estimate co-moving number densities corrected for 
the duty cycle ,^ri,c, i-e. ^ri,c = 7^ $mod, where 4™ indi- 
cates the cosmic time elapsed within each redshift bin, and 
Tjet indicates the lifetime of the radio jet. We find a global 
energy release at redshifts z= 1 — 3 (i.e., near the "Quasar 
Era"), dpE,kin/dt 



dpEMn/dt = 9 X lO^'^Ty ^ dEkineo/dteo erg s 



(1) 



where ry is given in lO'' yrs, and dE kin/ dtp^f^ is g iven in 
10*^° erg. For the redshift evolution of Willott et all (|2001h 
this corresponds to about 90% of the energy release at all 
redshifts, emitted during ~ 3.6 Gyrs. 

7.3. Overall feedback efficiency 

To quantify the overall efficiency of AGN feedback as ob- 
served in powerful radio galaxies, we will now compare the 
observed kinetic energy injection of the ionized gas within 
an activity phase with the possible total energy release 
through accretion onto the supermassive black hole. This 
simply corresponds to the relationship between the total 
energy equivalent of the rest mass of a black hole of a given 
mass and the kinetic energy of the outflows. 



77 = 5.6 X 10"^ 



EkinfiO tl Mg ^ 



(2) 
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Ekinfio indicates the energy injection in multiples of 10^° 
erg, ^7 the phase of AGN activity in multiples of 10^ yrs, and 
Mbh,9 the mass of the supermassive black hole in multiples 
of 10^ Mq. The efficiency i] is given in percent. 

Black hole masses for HzRGs are not easily derived 
directly. For an order-of-mass estimate, we will assume 
that the low-redshift relationship b etween black hole and 
bulge mass will approximately h old (*Mago rrian et al.lll998l : 
lFerraresel[2002l : iTremaine et~all[200 2: Hari ng k Rixll2004( l. 

The light profil es of HzRGs at z ~ 2 are regular and 

likely elliptical (|van Breugel et all Il998t iPentericci et al.1 
l200l| ). suggesting that these are bulge- dominated sys- 
tems- Adopting the stellar mass estimates of lSeymour et al.l 
(|2007t) for HzRGs, logMb„;o e 11 - 11.5, and the relation- 
ship of iHaring fc Rixl (|2004h . black hole masses are about 
log Mb?, ~ 8.2 - 8.8. This imphes that of order 0.1% of the 
energy equivalent of the black holes in HzRGs are released 
in kinetic energy of the outflows. This is a purely empirical 
estimate, and does not rely on any assumption regarding 
the efficiency of accretion or feedback. 

To investigate this question more globally, we can 
also relate the global energy density released by power- 
ful radio sources derived in pobs ~ 9 x 10^^ erg 
Mpc""^, wit h the local mas s dens ity of black holes, es- 
timated by lYu fc Tremai^ ()2OO20 . pbh = 2.9 ± 0.5 x 
10^ /ifo Mq Mpc~^, corresponding to a total energy equiv- 
alent of Pe,ki7i ~ 5.2 X 10^^ erg Mpc^^. This will serve as a 
strict lower limit on the efficiency with which the rest-frame 
energy equivalent of the black hole mass has been trans- 
lated into outflow energy in the past. We find an efficiency 
of about 0.2% with which the rest mass of the black hole is 
being translated into outflow energy, in agreement with the 
above estimate based on the masses of supermassive black 
holes in HzRGs, and with recent estimates comparing the 
optical and radio luminosity function of radio-loud AGN, 
but without th e direct observations of AGN-driven winds 
(jShankar et al.|[2008! ). 

The efficiency of AGN feedback in terms of the 
black hole accretion is an important parameter in 
cosmological models of galaxy evolution and struc- 
ture formation, which typically adopt a factor of_~ 
0.5% (e.g. [Sc annapicco & Oh 2004; Spr ingel et all l2005t 
iDi Matteo eTall 120051: [Chakrabarti et aP 120071 ). Given the 
uncertainties of factors of a few, and the fact that our esti- 
mates are based on the ionized gas only, the agreement is 
remarkable. This underlines the potentially significant role 
of jet-driven outflows for galaxy evolution. 



8. Summary and Conclusions 

We studied the kinematics and physical characteristics of 
the extended emission line regions in 3 powerful radio galax- 
ies at z~ 2, using integral-fleld spectroscopy of the rest- 
frame optical emission lines. Our results are consistent with 
what is expected from giant, jet-driven AGN outflows, sim- 
ilar to what we foun d in an earlier study o f the z~ 2.2 
galaxy MRC1138-262 (|Nesvadba et al.ll2006bl ). Analysis of 
the emission and continuum line morphologies, velocity 
maps, and flow orientation relative to the orientation of 
the radio jet does not favor alternative interpretations, e.g., 
merging, rotation, or gas infall. Such outflows may be com- 
mon in powerful high-redshift radio galaxies (HzRGs) dur- 
ing the "Quasar Era" . Overall, we find energetic outflows 



in all 4 2 powerful radio galaxies with extended jets 
for which we obtained rest-frame optical integral-field spec- 
troscopy. Comparison with previous longslit spectroscopy 
(jBaum fc McCarthy! |20000 suggests similar properties in 
> 10 HzRGs at z> 1. 

The emission line regions extend over 10 kpcxSO kpc, 
and are elongated along the axis of the radio jets. Most of 
the complex morphologies of these galaxies, previously ob- 
served with broadband imaging, appears to be due to line 
emission. The line-free continuum emission extracted from 
the data cube is relatively faint and compact (< 10 kpc). 
Assuming simple case-B recombination, we find of order 
10^° M0 of ionized gas, similar to the typical cold molecular 
gas content of HzRGs traced by CO line emission. From ob- 
served line ratios, we can constrain the physical properties 
of the gas. We detect [OHIJA4363 in a z~ 2 galaxy, allowing 
for a direct estimate of the electron temperature of ^ 10"* 
K. Electron densities are derived from the [SII]AA6716,6731 
line ratio and are ^ 500 cm^'^. From the measured Ha/H/3 
line ratio we estimate extinction of A(H/3)~ 1 — 4 mag. 

The velocity fields are very similar in all 3 sources 
and are reminiscent of two bubbles expanding back-to-back 
from the AGN. Velocity offsets near the center of the host 
galaxy are large, 700—1000 km s~^, but velocities within 
each bubble are relatively uniform. Line widths are of or- 
der FWHM~ 1000 km s^^, indicating strong turbulence. 
Larger widths in the central regions are likely due to over- 
laps between the bubbles. Comparison with the orientation 
of the radio jets suggest that this gas is outfiowing. 

With kinetic energies of order 10^" erg required to power 
the observed emission line kinematics, these outfiows may 
be sufficient to remove most of the interstellar medium of 
the host galaxies in a nearly explosive event. Observed ki- 
netic energies correspond to 0.2% of the energy equiv- 
alent of the rest mass of the supermassive black holes, in 
good agreement with AGN feedback prescriptions used in 
galaxy evolution models 

The following observations lead us to argue that the out- 
flows are dominated by the radio jets: (1) Geometry: Good 
alignment with the radio jets and elongation along the jet 
axis. The gas is blueshifted on the side of the approaching 
radio jet. (2) Size: The sizes of the emission line re g ions do 
not reach beyond the radio lobes. 'Nes vadba et al.l llOOT^) 
do not flnd extended outflows in galaxies with compact ra- 
dio sources observed with SINFONI. Emission line gas out- 
side t he radio lobes observed in Lya is kinematically quies- 
cent (jVillar-Martm et al.ll2003l ) . (3) Timescales: Dynamical 
timescales of the outflows are few xlO^ yrs, similar to typ- 
ical ages of radio jets. (4) Energies: Kinetic energies of the 
outflows correspond to about 10% of the jet kinetic energy. 

Comparin g with the z~ 2 radio luminosity function of 
IWillott et all (|2001h . we find that AGN winds fike those 
observed may release about 9x10^® erg s"^ Mpc^"^ in 
energy density from the massive host galaxies of power- 
ful radio sources at z~ 2. This energy may have a long- 
term impact on the subsequent evolution of the AGN host 
galaxy as well as nearby galaxies and help maintaining the 
"anti-hierarchical" properties of massive galaxies in spite 
of possibl e subsequent merging w i th companion galaxies. 
Following iNath fc Rovchowdhurvl (|2002D : iMcCarthv et al.l 
(|20'08l ). much of this energy may contribute to increasing 
the entropy and pressure in the environment of massive 
early type galaxies and clusters, and enhance the efficiency 
of dynamical interactions with galaxies in the environment 
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of HzRGs, such as ram pressure stri pping. Comparing with 
studies of ram-pressure stripping (jBoselh et al.l I2008D we 
find that the sateUite galaxies may aheady be red, poor in 
gas and rich in metals when they are being accreted onto 
the central galaxy. In spite of subsequent growth by up to 
~ 50% in stellar mass, we find that minor "dry" mergers 
with ram-pressure stripped satellites would not strongly in- 
fluence the observed evolutionary properties of the central 
galaxy. 
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Nesvadba et al.: Jet-driven AGN winds in Radio Galaxies in the "Quasar Era" 
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Fig. 8. top: [OIIIJA5007 emission line morphology of MRC0316-257 at z= 3.13 with contours showing the 1.4 GHz 
morphology, bottom left: Velocity map. Color bar shows the relative velocities in km s""'^, contours indicate the 1.4 GHz 
morphology, bottom right: Maps of the line widths, color bars show the FWHM in km a~^, contours indicate the 1.4 GHz 
morphologies. North is up, east to the left in all images. 
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Fig. 9. top left: [OIII]A5007 emission line morphology of MRC0406-244 at z= 2.42 with contours showing the 1.4 GHz 
morphology, top right: [OIIIJA5007 emission line morphology with contours showing the line- free rest-frame optical con- 
tinuum, bottom left: Velocity map. Color bar shows the relative velocities in km s"^, contours indicate the 1.4 GHz 
morphology, bottom right: Maps of the line widths, color bars show the FWHM in km s~^, contours indicate the 1.4 GHz 
morphologies. North is up, east to the left in all images. 
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Fig. 10. top left: [OIII]A5007 emission line morphology of TXS0828+193 at z= 2.57 with contoins showing the 1.4 GHz 
morphology, bottom left: Velocity map of TXS0828+193. Color bar shows the relative velocities in km s^^, contours 
indicate the 1.4 GHz morphology, top right: [OIII]A5007 emission line morphology with contours showing the line-free, 
rest-frame optical continuum, bottom right: Maps of the line widths, color bars show the FWHM in km s"^, contours 
indicate the 1.4 GHz morphologies. North is up, east to the left in all images. 
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Table 1. Sample and observation log 



Source 


z 


filter 


ToT 


seeing 


(1) 


(2) 


(3) 


(4) 


(5) 


0316-257 


3.13 


K 


100 


0.7"x0.4" 


0406-244 


2.44 


H 


90 


0.7"x0.6" 


0406-244 




K 


100 


0.8" X 0.6" 


0828-M93 


2.57 


H 


135 


0.7"x0.5" 


0828+193 




K 


170 


0.7"x0.5" 



Column (1) - Source ID. Column (2) - redshift. Column (3) - Total exposure time (on target) in min. Column (4) - Seeing. 
For TXS0828+193, the seeing was determined from a star within the field of view. Otherwise, we estimated the seeing from the 
standard star. 



Table 2. Emission Lines in MRC0316-257 



Line 


Ao 


Aobs 


z 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


[OIIIJ 


5007 


20665.4 ± 4.2 


3.1273± 0.0008 


1465± 90 


2.2±0.1 


[OIII] 


4959 


20459.6 ± 4.3 


3.1257± 0.0009 


1357± 94 


0.9±0.1 


H/3 


4861 


20072.2 ± 8.1 


3.1292± 0.0017 


902± 266 


1.2±0.4 



Column (1) - line. Column (2) - Rest-frame wavelength in A. Column (3) - Observed wavelength in A. Column (4) - Redshift. 
Column (5) - Rest-frame FWHM, corrected for instrumental resolution in km s~^. Column (6) - Line flux in units of 10~^^ erg 
s~^ cm~^. 



Table 3. Emission Lines in MRC0406-244 



Line 


Ao 




z 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




4861 


1.66532 ±0.00034 


2.42587±0. 00071 


1250±83 


1.9 ± 0.1 




OIII] 


4959 


1.69860 ±0.00034 


2.42528±0. 00069 


1524±93 


2.2 ± 0.1 




OIII] 


5007 


1.71495 ±0.00034 


2.42510±0. 00068 


1349±81 


6.6 ± 0.4 




Oil 


6300 


2.15777± 0.00045 


2.42503±0.00072 


826 ± 65 


0.74 ± 0.07 


Ha [NIIl 


6563,6583 


2.24879± 0.00045 


2.42647±0. 00060 


1756± 106 


9.75 ± 0.59 




SII] 


6716,6731 


2.30292± 0.00048 


2.42492±0. 00071 


1590± 102 


2.40± 0.16 



Column (1) - line. Column (2) - Rest-frame wavelength in A. Column (3) - Observed wavelength in A. Column (4) - Redshift. 
Column (5) - Rest-frame FWHM, corrected for instrumental resolution in km s~^. Column (6) - Line flux in units of 10~^^ erg 
s~^ cm~^. For Hq,[NII] AA6548,6583 and [NII]AA6716,6731 we give the values for the blended lines, respectively. 



Table 4. Emission Lines in TXS0828+193 



Line 


Ao 




z 


FWHM 


flux 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


H/3 


4861 


1.73883±0.00039 


2.57710±0.0008 


1448±112 


0.6±0.1 


[OIII] 


4959 


1.77422±0. 00036 


2.57778±0.0007 


1239±76 


2.1±0.1 


[OIII] 


5007 


1.79112±0.00036 


2.57723±0.0007 


1185±71 


5.9±0.4 


[01] 


6300 


2.2547±0.0006 


2.57821±0.0009 


1320±132 


0.5±0.1 


Hq + [Nil] 


6563,6583 


2.34705±0.00047 


2.57619±0.0007 


1176±47 


5.1±0.3 



Column (1) - line. Column (2) - Rest-frame wavelength in A. Column (3) - Observed wavelength in A. Column (4) - Redshift. 
Column (5) - Rest-frame FWHM, corrected for instrumental resolution in km s~^. Column (6) - Line flux in units of 10~^^ erg 
s"^ cm"^. For Hq,[NII] AA6548,6583 and [NII]AA6716,6731 we give the values for the blended lines, respectively. 



Table 5. Radio Properties of HzRGs 



Source 


Dang 


Dp/j,ys 


-Co. 1-1 


Tl.4GHz 


(1) 


(2) 


(3) 


(4) 


(5) 


MRC0316-257 


7.6 


58 


45.0 


0.94 


MRC0406-244 


7.3 


59 


45.2 


0.96 


TXS0828+193 


12.8 


103 


45.4 


1.7 


MRC 1138-262 


11.1 


92 


45.2 


1. 



For comparison we add MRC1138-262 from lNesvadba et al.l l|2006bl ). Column (1) - Source name. Column (2) - Maximal projected 
angular size of the radio source. Column (3) - Maximal projected size in kpc. Column (4) - Integrated luminosity between 0.1 and 
1.0 GHz in the rest-frame in erg s~^ cm~^, corresponding to ~ 0.001 — O.lx the kinetic luminosi ty. - Column (5) - Age of the 
radio jet for a constant expansion speed of O.Olc in 10^ yrs. Angular siz e distances are taken from lDe Breuck et al.l IpOOli '). radio 
luminosities are based on the VLA measurements of lCarilli et al.l l|l997f ). 
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Table 6. Outflow Properties of HzRGs 



Source Hopt Av Topt Ekiii.obs Ekin,ext E shell 

Sn (2) (3) (4) (5) (6) (7) 

MRC0316-257 21x9 670 3.1 44.4 

MRC0406-244 30x9 960 2.9 45.3 45.6 45.0 

TXS0828+193 28x10 800 4.1 45.3 45.3 45.0 



MRC1138-262 53x32 800 1 44.2 45 6 

For comparison we add MRC1138-262 from lNesvadba et a l' ('2006b'rc"olumn (1) - Source name. Golumn (2) - Extent of the optical 
line emission in kpc along the major and minor axis, respectively. Column (3) - Velocity offset between the blue and redshifted 
bubble in km s~^ in the rest-frame. Column (4) - Age of the outflow for constant Av and neglecting projection effects in 10^ yrs. 
Column (5) - Energy injection necessary to power the outflow in erg s~^ ('method 1', see S 15.51 for details), and mass estimates 
based on the measured Ha fluxes, uncorrected for extinction. This serves as a strict lower limit to the kinetic energy injection. 
Column (6) - The same as Column (5), but using extinction-corrected Ha fluxes. Column (7) - Energy injection necessary to 
inflate a hot bubble in erg s"'^, ('method 2', see i]5.5l for details). 



Table 7. Ha Luminosities and Ionized Gas Masses 



Source 




^ Ha 


£(Ha) 


£(Ha)"°'''' 


log(M/fc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0406-244 


-14.1 


-13.4 


44.5 


45.2 


10.6 


08284-193 


-14.2 


-13.4 


44.4 


45.3 


10.7 



Column (1) - Source ID. Column (2) - Measured integrated Ha flux in erg s~^ cm~^. Column (3) - Extinction corrected Ha flux 
in erg s"'^ cm~^. Column (4) - Integrated measured Ha luminosity in erg s~^. Column (5) - Integrated extinction corrected Ha 
luminosity in erg s~^. Column (6) - Ionized gas mass estimate. 



